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Abstract 

We argue that the statistical weight function W[p\ appearing in the 
McLerran-Venugopalan model of a large nucleus is intrinsically non-Gaussian, 
even if we neglect quantum corrections. Based on the picture where the nu- 
cleus of radius R consists of a collection of color-neutral nucleons, each of 
radius a <^ R, we show that to leading order in as and a/R only the Gaus- 
sian part of W[p] enters into the final expression for the gluon number density. 
Thus, the existing results in the literature which assume a Gaussian weight 
remain valid. 
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I. INTRODUCTION 



With the first colhsions at the Brookhaven Relativistic Heavy Ion Colhder, the experi- 
mental frontier in both energy and density has advanced. One of the goals of this machine 
is to search for the so-called quark-gluon plasma (QGP), a deconfined state of QCD where 
quarks and gluons form a sea of thermalized particles occupying a volume of space signif- 
icantly larger than that of a typical hadron. Theoretically, whether and when this occurs 
depends on the initial conditions of the system. That is, what is the distribution of the 
quarks and gluons immediately before the nuclei collide? 

One particularly fruitful approach to the determination of parton distribution functions 
in nuclei has been the McLerran-Venugopalan (MV) model The key observation 

underpinning the MV model is that for large enough nuclei and for small enough values of 
the longitudinal momentum fraction Xp, a new hard scale k^, corresponding to the large color 
charge per unit transverse area, enters into the problem. Consequently, when 3> Aqqj-,, 
it is argued that a classical treatment ought to provide a reasonable approximation to the 
gluon distribution. Three interrelated factors contribute to making the small-x^^ region 
conducive to a classical description. First, a large number of color charges contribute to the 
source, implying that the vector potential takes on large values within the nucleus. Second, 
the large vector potential corresponds to the existence of a large number of gluons, inviting 
us to apply the Weiszacker- Williams technique. And, third, the new large scale appears 
in the running strong coupling constant, q;s(k^), which is small for sufficiently large k^. In 
the classical treatment advocated in Refs. quantum mechanical expectation values 

are replaced by averages over a suitably chosen ensemble of color sources whose statistical 
weight function, W[p], was argued to take on a Gaussian form. 

Quantum corrections to the MV model were first computed in Ref . : they were found 
to be large due to contributions of order asln{l/Xp). The subsequent re-examination of 
the foundations of the MV model performed in Ref. led to the formulation of a set of 
renormalization group equations (RGE) describing the evolution of the weight function W[p] 
as one changes the separation scale between hard and soft partons @-14||. Early on in the 
development of the RGE it was conjectured that the MV model could be viewed as 
an effective theory at small x^, with the effects of the quantum corrections absorbed into 
a renormalization of the weight function Vr[p]. This conjecture has recently been proven 



by the calculation in Ref. |jT2|. While explicit solutions to the evolution equations for W[p] 
have not yet been obtained, it is known that a Gaussian distribution is not a solution to 
these equations [0. 

While these developments were taking place, we set out to address a different issue, 
namely the poor infrared behavior of the correlation functions in the MV model. The two- 
point vector potential obtained in Ref. grows like (x'^Y at large distances {x ^ Aq^j-,), 
signaling the onset of non-perturbative effects associated with confinement. In Ref. |TB[ we 
observed that since individual nucleons do not exhibit a net non-zero color charge, there 
should not be any long-range Aqqj-,) correlations between quarks. This requirement of 
color neutrality was cast into the form of a mathematical constraint on the two-point charge 
density correlation function. We found that the infrared divergence appearing in the MV 



model is completely absent when the color neutrality condition is enforced 



With the infrared divergences under control, it became clear that, at least in the classical 
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treatment of Ref. |T5|, the Xp dependence of the gluon distribution remained trivial, with 
XpdN/dXp independent of Xp. There are two fairly obvious possible sources of non-trivial x^ 
dependence, depending on the value of x^ under consideration. First, we can imagine taking 
Xp to be somewhat larger than the values specified in the MV model. In this situation, the 
longitudinal resolution of the gluons becomes good enough to start to probe the longitudinal 
structure of the Lorentz-contracted nucleus which they see. This avenue of investigation was 
pursued in Ref. |TB[ , where we developed a fully three-dimensional treatment of the classical 
gluon field of a large nucleus. Our view of the nucleus (radius R) as containing color-neutral 
nucleons (each of radius a) played a vital role in the calculation. Taking advantage of the 
smallness of the ratio a/R for a large nucleus, and assuming that the nucleus possesses 
a spherically symmetric distribution of color charge in its rest frame, we found that the 
additional Xp dependence manifests itself by the appearance of the combination q^+(Xj^m)^ 
in some (but not all) of the functions parametrizing the final result (see Eqs. (5.18)-(5.20) 
of Ref. [0). Here q is the transverse momentum of the gluon and m is the nucleon mass. 
Over most of the region where our extended treatment is valid, [XpTnY <^ implying that 
these corrections are small. 

The second possible source of non-trivial Xp dependence, namely quantum corrections 
proportional to asln(l/x^), was proposed even before the first calculation of these 

corrections in Ref. verified that they do indeed contain the necessary logarithms. The 
idea is that a resummation of such powers to all orders would convert the l/Xp appearing in 
the leading-order expression for dN/dXp into 1/x for some C. One of the goals of the 
studies conducted in Refs. [0-0 may thus be phrased as the determination of C. A priori it 
would appear that in order to determine the gluon distribution in the interesting nonlinear 
regime, we would require a complete solution to the evolution equations for W[p\ developed 
in Refs. And, furthermore, since a Gaussian is not a solution to these equations, it 



would seem that the general form of the gluon distribution obtained in Refs. |]5| JT5|JT6[| could 
be fundamentally altered once the non-Gaussian correlations are taken into account. Such 
an outcome would be at odds with related calculations based on different techniques 



which already possess general agreement with the MV model. The purpose of this paper 
is to investigate this issue. Once again it will prove valuable to incorporate the physics 
associated with confinement {i.e. the color neutrality of the nucleons) into the discussion. 
We will demonstrate that to leading order in the strong coupling and a/i? it is sufficient to 
consider only the (renormalized) Gaussian contributions to the gluon number density. All 
contributions from the non- Gaussian portion of W[p\ are suppressed by additional factors of 
Us and/or a/R. Thus, the MV model remains in general agreement with the various other 
approaches to small-x^^ physics referred to above, even if we employ a non-Gaussian weight 
function. 

In light of the arguments in favor of a Gaussian form for W[p\ put forth in Ref. [IH, the 
reader may be tempted to think that our conclusion is a trivial consequence of the central 
limit theorem. However, as we will demonstrate, the central limit theorem does not apply 
to this case. Instead, the additional factors which suppress the non-Gaussian terms relative 
to their Gaussian counterparts come about because of the fact that the individual nucleons 
are color neutral objects which do not have significant color-charge correlations at distances 
much greater than a. 

The remainder of this paper is organized into two main sections. In the first, we begin 
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with a review of the central hmit theorem to estabhsh the notation used in the subsequent 
discussion. We will pay particular attention to the conditions under which the central limit 
theorem is valid, and then show that for a large nucleus these conditions are violated. The 
essential observation is that any non-trivial longitudinal correlations, no matter how short 
range, spoil the Gaussian form of the weight function in the very large nucleus [A^^^ oo) 
limit. In the second half of this paper we will return to the three-dimensional framework 
for the MV model developed in Ref. [TO and consider the additional terms which would be 



generated in the presence of non-Gaussian contributions to W^[p]. Assuming only that the 
nucleus consists of color-neutral nucleons of radius a which possess no nontrivial correlations 
at separations much larger than a, we will demonstrate that no contributions beyond those 
already considered in Refs. |l5|, [T5|Jl6| are generated to leading order in and a/R. Thus, 



for a large enough nucleus, it is sufficient to know only the two-point function associated 
with even when the quantum corrections are taken into account. We close with a few 

words about the implications of our result. 

II. GAUSSIAN AND NON-GAUSSIAN DISTRIBUTIONS 

A Gaussian (or Normal) distribution in a random variable t usually comes about because 
of the central limit theorem. In this section we review the conditions under which this 
theorem is valid, and argue that all of these conditions are not, in general, satisfied by the 
color-charge distribution of a large nucleus, even when we omit radiative corrections. Hence, 
the statistical weight function W[p\ appearing in the MV model is generically not expected 
to be a Gaussian, even classically. 

A. Generating Functions 

Let T be a /^-dimensional random variable, with normalized distribution -P(t). For what 
follows it is more convenient to specify it by its moments: 

(TiT, ■■■Tk)= J d^'r TiT, ■ ■ ■ Tfc P(t). (2.1) 

These two descriptions are equivalent because the Fourier transform of P{t), namely 

P(A) = |rf^re'^^"P(T), (2.2) 
has a Taylor series expansion in terms of the moments: 

oo (_j)n 

pw = i + T.^—i-{ir-^r)- (2.3) 



n=l 



The moments defined by Eq. ( p.l| ) are reducible in the following sense. Suppose that the 
distribution P(t) factorizes as 



D 



P{r) = n Pin), (2.4) 



i=l 
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where p{Ti) is some single variable distribution, i.e. suppose that the components of the 
vector T are completely unrelated to each other. In this case, the moment {TiTj ■ ■ • r^) 
factorizes into {Ti){Tj) ■ ■ ■ (r^). More generally, however, the components of r will not be 
completely independent of each other. Then, the factorization just described will fail: for 
example, (TjTj) ^ The difference between these two quantities provides a measure 

of the amount of interdependence between the components of r. This information may be 
conveniently organized in terms of the (irreducible) cluster moments {{TiTj ■ ■ ■ Tk)), which are 
recursively defined through the relations 

in) = {{n)), (2.5a) 
(r.r,) = ((r,r,)) + ((r,))((r,)), (2.5b) 
{r^r,n) = ((r,r,r,)) + ((r,))((r,r,)) + ((r,)) ((r,r,)) + ((r^)) ((r,r,)) + ((r,)) ((r,-)) ((r,)), (2.5c) 

{TiTjTkTi) = {{TiTjTkTi)) 

+ {{ri)){{Wkn)) + {{r,)){{rirkri)) + {{n)) {{riTjn)) + {{n)) {{nrjTk)) 

+{{nrj)){{nTi)) + {{nn)){{r,n)) + {{nn)) {{r^n)) 
Hn)){{rj)){{nn)) + {{r,)) {{n)) {{r,n)) + {{n)){{n)){{r,n)) 
Hr,)){{n)){{nn)) + {{r,)){{n)){{nn)) + {{n)) {{n)) {{nr,)) 

Hr^)){{r,)){{r,)){{n)). (2.5d) 

The generalization to higher orders is obvious. It turns out that the cluster moments are 
generated by InP(A), in the same way that the ordinary moments are generated by -P(A): 

oo (_j)n 

lnP(A) = E^(((T-A)")). (2.6) 

n=l 

A simple example is provided by the following Gaussian distribution: 

All of the odd moments (reducible and irreducible) of Eq. ( p.7| ) vanish. Of the irreducible 
moments, only the two-point function {{TiTj)) is nonvanishing. All other cluster moments 
are zero, and the even reducible moments are expressible entirely in terms {{tiTj)). 



B. Central Limit Theorem 



Suppose the random variable r is independently sampled N times, yielding the values 
T-ip) (p = 1,2,..., N). The central limit theorem asserts that when iV 3> 1, the mean value 
of these measurements. 



obeys a Gaussian distribution with 

m = {{n)) = in) 

and 



{{titj)) 



N 



In other words, the distribution of t is given by 



W{t) 



detn 
W 



where 



-1 



AT 



(2.9a) 



(2.9b) 



(2.10) 



(2.11) 



The proof of the central hmit theorem rehes crucially on the independence of samplings: 
the joint distribution Q of the sampled variables r*-^-* must be factorizable, i.e. 



Q{r 



(1) T-(2) _ 



(7V)^ 



N 

np(r(^)). 

p=l 



(2.12) 



In terms of Q, the probability W{t) for finding the mean value t is given by 

W{t) = J d'^T^^') J d''T^''> ■■■ j d'^T^''^ Q{t^^\t^^\ . . . , r(^)) 5^ (t - E T(f)/iv') . (2.13) 

p=i 

Using the distribution given in ( p.l2| ), the moment generating function is simply 



W{\) 



d^t e-'^ 'W{t) 

N 



i\ t ^ 



This leads to the cluster moment generating function 

InW^(A) = N\nP{X/N) 

-zrnir-xr)) 



(2.14) 



E 

n=l 



N' 



n-l 



(2.15) 



where we have applied Eq. ( p.6| ) to obtain the second line. Now consider the ^ 1 
limit. If we approximate (|2.15| ) by just the first term, we obtain iy(A) = exp(— zA ■ (r)), 
corresponding to W{t) = 6^{t — (r)). If instead, we include the first two terms in the 
infinite series, we get 



W{X) = exp(-a-(T) - XiXj{{nrj))/2N 



(2.16) 



corresponding to the Gaussian distribution (|2.10|) . Deviations from this Gaussian are sup- 
pressed by the additional powers of 1/iV present in the higher order terms. 
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C. Non-Gaussian Distribution 



The proof of the central hmit theorem depends critically on the assumption of indepen- 
dent sampling. If (|2.12|) is violated, then the proof fails. Even a small correlation between 
successive samplings is sufficient to destroy the conclusion, as the following example illus- 
trates. 

Suppose we introduce a "nearest-neighbor" correlation by replacing (|2.12 ) with 



Q{r 



(1) _(2) _ 
,1 , . . . , I 



N-l f(^(m) _(m+l)N 
1+^ /IT ,T 



P(t('")) P(r(™+i)) 



N 



(2.17) 



To ensure that ( |2.17| ) is properly normalized, we further assume /(r'^™-', t*^™+^)) to be sym- 
metric in its arguments, and that 



The moment generating function for the distribution ( p.l3|) is now 



(2.18) 



W{\) = P{\/N) 



N 



1 + 



{N - 1)/(A) 



P{X/N) 



(2.19) 



where 



/(A) = I d^T^'^ I d^T(2) e-^^-l-^'^+^^'^l /(t(i),' 



.{2)^ 



(2.20) 



Note that /(O) = 0. Unlike Eq. ( |2.15| ), the irreducible moments generated by 



In W{X) = N\n P{X/N) + In h + 



(N - l)f{X/N) 



PiX/N) 



(2.21) 



do not approach a linear function of A as ^ oo. In fact, we have 



lnW{X) = -zA-(r) + In [1 - + 0{1/N), 



where 



d 



/(A) 



A=0 



(2.22) 



(2.23) 



For (p ^ 0, Eq. (|2.22|) is nonlinear in A. This distribution possesses non-zero values for all 
of its irreducible moments, independent of A^. Including the 0{1/N) contributions does not 
change this conclusion. The corresponding distribution W{t) is no longer a pure Gaussian: 
the central limit theorem does not hold in this case. While it is straightforward to determine 
the (new) form of the distribution function for this particular example, in general this task 
will not be so easy. Fortunately, all we need to know for what follows is the fact that the 
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resulting distribution is fundamentally non-Gaussian when the successive measurements are 
(even mildly) correlated. 

It is useful to see a bit more about how the correlators change when we relax the assump- 
tion of independent measurements. With independent sampling, as in Eq. (|2.12|) , correlators 
such as {r-^^Tj'^^ ■ ■ ■rjf'^) factorize into groups with the same superscripts. For example, 



= {n){T,){ruTi), (2.24) 



assuming that p, and s all take on distinct values. Of particular interest in connection 
with the MV model will be the special case 

{rl'\f^) = 6^\nT,) + (1 - 5^^)(r,)(r,). (2.25) 

In contrast, if the measurements are not independent, such as in ( |2.17| ), then the factor- 
izations obtained in Eqs. ( p. 241) and (|2.25| ) are no longer valid. Additional non-factorizable 
correlators involving different superscripts will, in general, be present, signaling the break- 
down of the central limit theorem. 



D. Color Charge Distribution in a Large Nucleus 

Let p"(a;|| , x) be the color charge density of a large nucleus with color a and atomic weight 
y4 1, at the transverse position x and the longitudinal position (a precise definition of 
is given in Eq. ( |3.3|) below). Let 

/oo 
dX||p"(X||,a;) (2.26) 
-oo 

be the accumulated two-dimensional charge density. On account of Lorentz contraction of 
the longitudinal size of the relativistic nucleus, this is expected to be the relevant quantity 
governing the gluon distribution, provided that only valence quarks are included in p, and 
provided that the longitudinal wavelength of the gluon is much bigger than the longitudinal 
size of the relativistic nucleus. These are two of the basic assumptions underlying the MV 
model g-i. 

The typical size of P2{x) for a nucleus is of order A^^^ times the typical size of P2{x) for 
a single nucleon. If we identify A^/'^ with the number of measurements of the previous 
subsections, and the accumulated charge density P2{x)/A^/^ with the mean value (with 
i —> {a,x}), then the central limit theorem would force W{t) (that is, iy[p2(a^)]) to be 
Gaussian in the limit of large A^^^, provided that the conditions needed to prove the theorem 
are obeyed. With this mapping, the superscript p in the random variable r/^^ corresponds 
to and r/^"* corresponds to p^{x^\, x). A summary of these connections appears in Table |. 

As discussed in Sec. Ill C| , a necessary condition for the central limit theorem to hold is 
that the correlators {rj^^^j ■ ■ ■t^^'') must factorize into groups with identical superscripts. 
In particular, consider the translation of Eq. ( p.25|) into the MV model. The color neutrality 



condition tells us that (p"(X||,a;)) = 0. Therefore, we are left with 
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{{p\x,,x)p\x\,x'))) = {p'^ix„x)p\x[,x')) 

= 5{x,-x[){pl{x)p\{x')). 



(2.27) 



This is effectively the form of correlator argued for in Ref. [|l|, and employed extensively in 
the literature ||^,|l3,|ll,|3,||-|8l . 

However, there are two reasons why it is desirable to go beyond the Gaussian approxi- 
mation represented by Eq. ( p.27| ). First, at very small x^ the quantum corrections in the 
MV model become large. If we wish to incorporate these corrections via the RGE analysis 
of Refs. 1 we must go beyond the Gaussian approximation, since a Gaussian is not a 



solution to these equations [0]. Second, if we go to somewhat larger values of x^, where the 
quantum corrections may not be too large, the presence of longitudinal correlations between 
quark charges within a nucleon and the fact that the longitudinal size of a nucleon is not 
exactly zero begin to have an effect, even classically. In Ref. |T^, we employ the form 



ip''{x\\,x)p\x'x'))) (X d^^l^ix^^ - x'^ 5^{x - x') - C{x\\ -x[]x 



X 



(2.28) 



where C{x\\ — x'^x — x') is a reasonably smooth function parameterizing the mutual cor- 
relations between pairs of quarks. This function enforces the color neutrality condition. 
The nonfactorizability of Eq. (|2.28|) tells us that for a nucleus with a non-zero longitudinal 
thickness we should not expect any single layer to be color neutral by itself. Equivalently, 
we may view the color neutrality condition as forcing us to have non-trivial longitudinal 

p!6[| , all of the higher-order cluster moments of the distribution corre- 



correlations. In Ref. 



sponding to Eq. (p.28|) are taken to vanish, implying that W[p"'{x\\ \ x)] is Gaussian. However, 
Eq. ( |2.27| ) is violated, and the corresponding 2-dimensional distribution iy[p2(ic)] is not a 
Gaussian, because the central limit theorem does not apply to this computation. Nev- 
ertheless, the result obtained in Ref. |T^ remarkably turns out to have exactly the same 
form as if iy[p2(^c)] had been Gaussian: for sufficiently small x^, it matches the results of 
Refs. [H,|l5l. In the next section we will see that this outcome is quite general, by considering 
the consequences of allowing for a non-Gaussian weight function. 



III. EFFECT OF NON-GAUSSIAN CONTRIBUTIONS 

We have just demonstrated that we do not expect the statistical weight function TV [^2(2?)] 
to have a Gaussian form, even in the purely classical case. On the other hand, a Gaussian has 
been frequently employed in the literature dealing with the MV model [|I]-|,|T^,|r^,^,^^-^ 



In this section we will argue that for a large enough nucleus, the results which would be 
obtained with a non-Gaussian weight function are identical to those which follow from a 
Gaussian weight function to leading order in «<, and A^/^. That is, the new contributions 
introduced by the non-Gaussian weight are suppressed by additional factors of as and/or do 
not contain as many factors of A^/^. Our conclusion is a consequence of the color neutrality 
of the nucleons: all non-trivial correlations are limited to length scales of order a ~ Aq^j-,. 
Beyond the color neutrality of individual nucleons, we will assume nothing about the form 
of Vl^[p(x)]. Our discussion is framed in terms of the three-dimensional extension of the MV 
model introduced in Ref. [ffBl, with the addition of non-zero higher-point cluster moments. 
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A. Diagrammatic Representation of the Gluon Number Density 



Imagine a very large nucleus as viewed in its rest frame. The current corresponding to 
this situation may be written in the simple form 

where we employ the subscript "r" to denote rest frame quantities. In this frame the Yang- 
Mills equations for the vector potential possess the "obvious" time-independent Coulomb 
solution. Furthermore, since only ^ 0, we have BqA^ = d ■ = the Coulomb solution 
is the same as the covariant gauge solution in this frame. Thus, we conclude that even 
when we boost to the lab frame, it is natural to begin with the covariant gauge solution 
for the vector potential |jl6|. Indeed, it has been observed that although the expression 
for the gluon number density is most easily written in terms of the vector potential in the 
light-cone gauge (see Eq. ( |3.10|) below), it is nevertheless easiest to work in terms of the 
covariant gauge expression for the current []T^|T3|,|T^^|. 



In the lab frame, where the nucleus is moving along the +z axis with speed /3, we take 



the source to be of the form [IC 



J+ = - p{x\^,x)] J- = ^J+- J = 0. (3.2) 



The longitudinal coordinate x,, is defined by 



e 2 

The parameter 



= - X — - x~^. (3.3) 



'2(1-/3) 



/3 



(3.4) 



measures how close the source is to being exactly on the light cone.Q As explained in Ref. \TE 
we prefer to work with natural (order unity) quantities, and keep track of all small and large 
parameters explicitly through e. At the end of the calculation we let e — 0. Interestingly, 
in terms of Xy, the function describing the source is still spherical: the Lorentz contraction 
that shrinks x~ towards zero in the lab frame is exactly compensated by the factor 1/e. 
Reflecting this fact, we will frequently use the notation x = {xf, x). 

As explained above, it is natural to first solve the Yang-Mills equations for the vector 
potential in the covariant gauge, using the covariant gauge source ( p.2| ), and to subsequently 
transform the result into the light-cone gauge = 0. The result of this procedure may be 
written in the form 



^We define the light-cone coordinates to be = —x^ = {x^ it x^)/\/2- The transverse coordi- 
nates x^ and form a two-vector which we write in bold-face: x. Our metric has the signature 
(— , -|-, -|-, -|-). Thus, the scalar product in light-cone coordinates reads g^x^ = —q~^x~ — q~x^ + q-x. 
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IQH J-oo J 

x(ri/ d%G{x,,-^,,;x-^,)Y[p{^^)p{^2)---p{U]], (3.5) 
where is the momentum conjugate to Xy, i.e. 

roc 

A^{q^f.,x)= / exp(— igiiXii) as). (3.6) 



The power counting rules we will present in Sec. [1IIB| are simplified by using A{qf, x) instead 

of A(X||, x). 

Although Eq. ( p. 51) appears complicated at first glance, it is easily understood in terms 



of the diagrammatic representation introduced in Ref. |]T6[. The light-cone gauge vector 
potential is a non-linear function of the source, containing all possible "powers" of p. Note 
that the series begins at order g, and that each additional occurrence of the source adds 
a factor of g^. Fig. ^ illustrates the first few diagrams corresponding to the series. The 
mth order diagram consists of m copies of the source, represented by the circled crosses. In 
Eq. ( p.5| ) they appear as the nested multiple commutator 



p(6)p(6)---p(e„ 



p(ei),p(6) ,p{C3 



(3.7) 



To each source we attach a propagator (Green's function) connecting the source point 
(^llfc,Cfc) to the point (X||^,a;): 



X 



life 



e,iJ^ + (^-6 



(3.8) 



Here x is the transverse position at which we wish to know the vector potential. The index 
labelling the uppermost propagator on the diagram represents the derivative indicated in 
Eq. ( |3.5|) . Finally, we have an ordered integration over the longitudinal variables Xy^: 



dx II 



dx II 



dxu 



(3.9) 



Mathematically, this integration results from the gauge transformation to the light-cone 
gauge from the covariant gauge. Physically this integration corresponds to the final state 
rescatterings which would be present in a computation of the gluon number density based 
entirely on the covariant gauge p8|-plj|. The ordered integration is represented by the dots 
on the vertical line. The dots are to slide up and down the entire length of the line without 
passing one another. Although we have written all of the integrations over an infinite range, 
in practice the source provides non-zero contributions only over a region of size R, the 
nuclear radius. 

The next step is to connect the gluon number density to the two-point correlation function 
for the vector potential [jl6|,^: 



dN 



1\\ 



dq«d'^q Att^ 



J d'x J d'x'e''^<^-^'\At{q,;x)A^{-q,;x')). 



(3.10) 
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Intuitively, the result in Eq. ( |3.1CI|) may be understood by envisioning the expansion of the 
vector potential in terms of creation and annihilation operators and recognizing that (AA) 
contains the number operator PJ^. The vector potential appearing in Eq. ( |3.10| ) is in the 
light-cone gauge. This choice reflects the fact that the intuitive picture of the parton model 
is most transparently realized in the light-cone gauge p9|-[5^. Strictly speaking, the angled 
brackets in Eq. (|3.1CI|) represent a quantum-mechanical expectation value. In the MV model, 
we make a classical approximation to this quantity by performing an ensemble average with 
an appropriate weighting function iy[p(x)]. 

Given a specific form for we may (in principle) evaluate Eq. ( |3.1(J| ) by inserting 

two copies of ( p.5[ ), performing the required average term-by-term, and summing the result- 
ing series. Diagrammatically, the quantity {AA) may be represented by drawing all possible 
pairs of diagrams for a single A, with the understanding that all possible contractions should 
be performed (see Figs. ^|-^ for the first three orders in this expansion). For a Gaussian 
weight function, only pairwise contractions appear. On the other hand, for a non-Gaussian 
distribution, contractions connecting three or more sources must also be considered. Be- 
cause the nucleons are neutral, (p(x)) = 0. Hence there must be no uncontracted sources 
left over in any diagram. 



B. Power Counting Rules 

To proceed further, we must make some reasonable assumptions about the correlation 
functions corresponding to the moments of Vr[p(x)]. Let us parameterize the m-point 
correlation function (m > 1) by 

{{p{Xi)p{x2) ■ ■ ■ p{Xm))) = Smi—y^Xj) V„,{xi-X2\ X2-X3; . . . ] X„_i-f„). (3.11) 

The exact choice made for the m—1 difference coordinates is arbitrary and unimportant to 
our argument. What does matter is the fact that the source corresponds to a large nucleus of 
radius R which is in turn composed of A nucleons, each of radius a. The nucleons themselves 
are color-neutral. Because of confinement, we do not expect the field to be correlated at 
distances much greater than a: what is happening inside one nucleon is largely independent 
of what is happening inside of the others. Thus, the function ought to be small unless 
for all pairs of points. Furthermore, the center-of-mass coordinate ought to 
point at a position somewhere inside the nucleus {i.e. it should have a magnitude < R) in 
order for to take on a non-negligible value. 

The above physical considerations are sufficient to allow us to determine the order of 
magnitude of an arbitrary diagram in terms of powers of a, and R/a {A^'^). The powers of 
the coupling are simple. Recall that each of the diagrams for the gluon number density ( |3.10| ) 
are formed by gluing together two copies of the expansion for the vector potential ( p. 51) . 
Since ( |3.5| ) contains only odd powers of the diagram representing the contribution to ( |3.10| ) 
with a total of j sources contains the factor g^^'"^ (or, equivalently Oi{~^). 

Now let us consider the integrations which go into the computation of the contributions 
to (yl^(g|j; a;)yl"(— gy; x')). The integrand for a given diagram will contain several propagators 
plus factors of Sm and depending upon how the sources are contracted. As described 
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above, Vm contains the length scale a whereas Sm contains the length scale R. The Green's 
functions contain no intrinsic scale of their own. What we would like to know is how many 
powers oi R/a ~ A^^^ are generated when we perform the required integrations, which range 
over all possible locations of the sources as well as the entire length of the two vertical lines. 
To first approximation, the function Sm in Eq. is essentially a constant when we 



stay well inside the large nucleus since the individual nucleons are identical insofar as the 
strong interaction is concerned. Thus, the quantity {A'^lqf, x)A'^{—qf x')) computed from 
the expansion given in Eq. ( p.5| ) is of the form 

{R/arF{q,; x - x') ^ {A'/'Y F{qf, x - x') (3.12) 

where F is a dimensionless function. For a given diagram, the power c is equal to the 
number of independent clusters formed when the vertical lines are removed. Two clusters 
are independent if they can be represented by a planar diagram (no crossings). For example, 
c = 1 for the diagrams in Figs. 0, |^, and ^(f)-^(k), whereas c = 2 for the diagrams in 
Figs. 1(a) J(e). 

The origin of the factor of A^^'^ for each independent cluster is the divergence which 
appears in some of the integrations when we take Sm = const.: these are the integrations 
which get cut off at the scale R instead of the scale a. To see this, imagine that we have 
changed the integration variables to a set of sum and difference coordinates. The sum 
(center-of-mass) coordinates describe the positions of the clusters (whether or not these 
clusters are independent).^ The difference coordinates describe cluster-cluster separations 
as well as the internal separations of the components of each cluster. All of the difference 
integrals are finite and so inherit the scale a associated with Dm- On the other hand, the 
center-of-mass integrations are insensitive to a: so long as the cluster is located somewhere 
inside the nucleus, the integrand can be significant. Hence, these integrations produce factors 
of R instead of a. In the absence of the ordering of the locations of the points on the vertical 
lines, then, we would obtain a factor {R/aY for k clusters. What the ordering does is to 
force clusters which cross each other in nonplanar diagrams to maintain a relative separation 
which does not greatly exceed a. Thus, if only c clusters are free to move independently of 
each other, the factor we obtain will be only [R/aY ~ [A^/'^Y- 

Finally, we note that every diagram will pick up an additional factor {R/aY ~ A'^^^ when 
the Fourier transform indicated in Eq. ( |3.1(J| ) is performed. The origin is the same as above: 
the integration over x—x' gets cut off at the scale a whereas the integration over {x+x')/2 
is cut off at the scale R. Assembling all of the pieces, we have the result that the diagram 
containing j sources and c independent clusters contains the factor ai~^{A^/'^Y^'^ ■ 

With our power-counting rules in hand, we may now classify the relative contributions 
of the diagrams in Figs. ^]-^. The lowest order diagram is given in Fig. 0, and contains 
two sources and a single cluster. Hence, its contribution is of order asA. The third-order 
diagrams of Fig. |^ have a total of three sources tied together as a single cluster. Their 



^Note that transverse position {x + x')/2 of the clusters is not integrated over when computing 
{A^{q\\]x)A'^[—q\\;x'Y- Hence, the only center-of-mass integrations present are longitudinal in 
nature. 
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contribution is therefore order a'^A and is suppressed by a power of as relative to the leading 
order diagram. Finally, let us look at the fourth-order diagrams of Fig. |[ According to our 
power-counting rules, diagrams (a)-(e) all contribute at the a^A'^^^ level (four sources and 
two independent clusters), whereas diagrams (f)-(k) produce only a^A. We imagine that 
we have a very large nucleus such that a^A^^^ is of order 1 or greater. Thus, we consider 
diagrams (a)-(e) to be as important as the leading order, whereas diagrams (f)-(k) are 
subleading, being suppressed by one factor of a/i? ~ A~^^^. 

The generalization to higher orders is obvious. Clearly, the leading diagrams are all 
ladders: they contain only two-point contractions which do not cross, thus producing the 
maximum number of independent clusters for a given number of sources. All other diagrams 
are suppressed: either they contain extra factors of as or they have fewer factors of A^^^. 
Therefore, the A and as dependence of the leading contributions to the final result must 
look like 

X, = CoasA(l + cia^^Vs + caa^A^/^ + ...). (3.13) 

But this is exactly the form obtained by expanding the expressions contained in Eqs. (5.18)- 
(5.20) of Ref. [0. In fact, the set of diagrams described above is precisely the set of diagrams 
considered in Ref. |T^. Moreover, it is also the set of diagrams whose two-dimensional 
reduction corresponds to those terms which were retained in Refs. p|,p!5[|. These contributions 



are easily resummed into a relatively simple "exponential" P JT5|JI6|] . Thus, we see that even 



if the weight function iy[p(x)] is non-Gaussian, we obtain the same result to leading order 
as if we had chosen to use a Gaussian weight function instead. 



IV. IMPLICATIONS 

We have just demonstrated that the non-Gaussian portions of the statistical weight 
function FT [p(a; )] in the MV model do not contribute to the gluon number density at leading 



order. On the other hand, the results of Ref. [12| suggest that the quantum corrections to 



the MV model may be incorporated by solving the RGE for the weight function and using 
this (non-Gaussian) result to perform the ensemble averages required to compute the gluon 
number density. A synthesis of these two conclusions has the following implication: it is 
sufficient to solve the RGE for the new value of the two point charge-density correlation 
function, and to use this renormalized correlator as defining an effective Gaussian weight 
function to be input to the MV model. The separation scale between hard and soft partons 



appearing in Ref. ||T2[ may be recast as a dependence of the renormalized correlator on 
Xp. In terms of Eq. ( p.28| ), this dependence can potentially show up in the detailed shape 
of the smooth function C{x — x'), as well as in the prefactors which were suppressed in 



writing down Eq. ( ^.28| ). In terms of the expansion ( |3.13D , this means that the Cj's would 
be x^-dependent. Since the quantum corrections presumably incorporate the physics of the 
Yukawa cloud believed to surround the nucleons, it is possible that the detailed relation 
between a and R could change somewhat, again in an x^-dependent way. 

In Ref. 1^ it was argued that not only should the RGE-improved MV model apply to 
nuclei at sufficiently small values of x^, but that it should also apply to hadrons as well, at 
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even smaller values of . Our argument does not apply to this situation. For a hadron we 
would effectively have a/R = 1, meaning that all of the suppressed contributions we have 
dropped are no longer unimportant. In this situation the full non-Gaussian solution to the 
RGE has relevance and must be taken into account. 
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FIGURES 




FIG. 1. Diagrammatic representation of the series expansion for the hght-cone gauge vector 
potential, Eq. (|3.5|). The circled crosses denote the positions at which the sources are being 
evaluated. The dots represent the ordered integrations coming from the gauge transformation into 
the light-cone gauge. A propagator connects each of these points to the sources. 
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FIG. 2. The lowest-order contribution to the gluon number density, Eq. ( ^.10 ). According to 
the power-counting rules described in the text, this diagram is proportional to OgA. 
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FIG. 3. The contributions to the gluon number density ( 3.1C| ) at 3rd order in p. According to 
the power-counting rules described in the text, both diagrams are proportional to a^A. 
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FIG. 4. The contributions to the gluon number density ( ^.lOf ) at 4th order in p. Diagrams 
(f)-(h) are non-planar, whereas diagrams (i)-(k) represent the (non-Gaussian) contribution from 
an irreducible 4-point charge density correlation function. According to the power-counting rules 
described in the text, diagrams (a)-(e) are proportional to a\A^/'^ whereas diagrams (f)-"(h) are 
proportional to a\A. 
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TABLES 



TABLE I. Mapping between the notations used in the discussion of the central hmit theorem 

and those associated with the MV model. 

central limit theorem MV model 

i {o-]x} 
N AV3 

W{t) W[p2{x)] 
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